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Abstract. The lowest-order induced soft gluon radiation processes of heavy 
quarks have been analyzed to quantify the dead cone effect. This effect is most 
likely expected to suppress significantly the energy loss of charm quarks passing 
an amorphous color-neutral deconfined medium, as has been concluded from 
recent experiments at RHIC 
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1. Introduction 

Radiative energy loss of quarks traversing strongly interacting matter is considered 
as a promising probe of the quark gluon plasma. Indeed, the recent heavy-ion col- 
lision experiments at RHIC show a distinct suppression of transverse momentum 
spectra of hadrons composed of light quarks [ U thus evidencing the transient cre- 
ation of a medium with large stopping power caused by a high gluon density. In 
contrast, open charm mesons [ |2] seem to suffer only a tiny, if any, energy loss 
[E]. As noted in [0j prior to this experimental finding, the gluon emission pat- 
tern of heavy quarks experiences a modification, known from jet physics as dead 
cone effect, which reduces the radiative energy loss. While large efforts are de- 
voted to the study of the energy loss of light quarks at asymptotic energies (cf. [ 
EHSH3|Hlini[ini[ni[I2[T3|and further references therein), the treatment of heavy 
quarks [O El is m progress. 

Similar to the fairly complete treatment of the radiative energy loss of light high- 
energy quarks we consider here the lowest-order induced gluon emission of a heavy 
quark. It is our goal to take exactly into account the finite mass, the finite energy 
and correct kinematics (i.e. not only soft gluon emission), the full color algebra 
and all spin/polarisation effects. In doing so we restrict ourselves to the single 
scattering and one-gluon emission process. The importance of the single scattering 
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(a) Mi, 0,0 
(pre emission) 



(b) Mi, 1,0 
(post emission) 



(c) Mi, 0,1 (three gluon) 



Fig. 1. Diagrams of single scattering with one-gluon emission classified by the 
notation scheme M„, m ,; of [Ej- 



process is emphasized by various aspects. Among of them data from RHIC, which 
suggest thin plasmas, imply that only a few scatterings are taking place; single 
scattering is the first, and according to [|7| the most important. The existence of an 
analogue to the Landau-Pomeranchuk-Migdal (LPM) effect in QCD, which reduces 
the effective number of scatterings further, is another indication of the significance 
of this process. Moreover the fast convergence of the opacity expansion [0, where 
the first order describes the radiative energy loss considerably well, suggests the 
relevance of this basic energy loss process. The discussion of the QCD analog of the 
Ter-Mikaelian (TM) effect [E] is also based on the lowest order in opacity [ II Tj 1 . 

The article is organized as follows. In Section 2 we describe the potential model 
and the concept of radiation amplitude and dead cone factor. Numerical results are 
then discussed in Section 3. The validity of the potential model is commented on 
in Section 4. The last part summarizes conclusions. 

2. Single Scattering in the Potential Model 

Following [ we consider the processes depicted in Fig. 1. The crosses stand 
for external potentials which, according to [ 13,, represent the medium. A color 
charge, here a quark, moving through a background field of quarks and gluons is 
mimicked by scattering on target particles modeled by static screened potentials 

localized at the space points Xi, V^ B (q) — gT% B — r, q is the momentum 

q + fi 

transfer and /i is the Debye screening mass; the coupling g is not varying here. 
The term T\ B = ^ B C^'''^ m \A stands for the original color structure of the target, 
where T a are generators of SU(3) co i or , \ the quark color states and A,B = 1 ... 3 
are color indices. In order to allow for the neglect of radiation contributions from 
targets in the realistic situation where the potential in Fig. 1 is substituted by 
a dynamical target, and thus for a realistic model, one requires a particular gauge 

lr The results of f !17| confirm the conjecture f llril that the asymptotic gluon mass, which is also 
important for describing phcnomcnologically the quark-gluon plasma f 1181 . determines the TM 
effect. 
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choice, the A + -gauge. The gauge choice is given by the condition A + = A° + A z = 
in light-cone coordinates. 

In the soft radiation limit the conditional probability of one-gluon emission 
under the condition that an elastic scattering has occurred is given by the ratio of 
the inelastic and elastic differential cross sections 



dp _ dgjnd _ |M incl | 2 d 3 fc 
~ d^ " \Mdf 2 W (2^) 3 ' 

where the overline indicates an average over initial polarization and color configura- 
tions and a sum over all final polarizations and colors. Thus the sufficient quantity 
to compare different radiation patterns in such scattering processes is the radiation 
amplitude |i?| 2 = |Afi no i| 2 /|Af c i| 2 . If the momentum of the gluon is small, compared 
to the momenta of the other particles, a factorization into a part describing the 
emitted radiation and an underlying part for the elastic scattering can be carried 
out on the level of the matrix elements. 

The probability, i.e. the radiation amplitude, for gluon emission off heavy quarks 
in projectile direction is strongly reduced at radiation angles 9 < 9q by a dead cone 
suppression factor according to F 2 [0] with 

which corrects the matrix elements derived in the approach of Gunion and Bertsch 
[ 19 for massless, high energetic quarks. Below 9o = m/E, where E means the 
energy of the incident heavy projectile quark, radiation becomes suppressed due to 
the mass m of the projectile, which introduces a new scale relevant for the radiation 
pattern. However, the suppression factor F is restricted to small angles. Moreover, 
this estimate is valid for soft gluons, that means it focuses on abelian diagrams only, 
and cannot properly handle intermediate gluon energies. These shortcomings can 
be overcome only by a numerical computation, which is the main motivation for 
the present work. 

3. Numerical Results 

The graphical elements depicted in Fig. I correspond to objects in our CH — h code 
[ 120). Therefore, the matrix elements are evaluated numerically. All interference 
effects are exactly dealt with when squaring sums of matrix elements and perform- 
ing averaging/summation over spins/polarisations. Since spinor calculations do not 
reveal significant modifications we present results of the scalar QCD case and com- 
ment later on spin effects. 

3.1. Abelian Region 

To highlight the role of the non-abelian diagram lc with the triple gluon vertex we 
consider first the "abelian part of QCD". In Fig.|2the radiation amplitude obtained 
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Fig. 2. Radiation amplitude |i?i| 2 as a function of 9 for m < \q±\ (light projectile 
quark, ma = 0.007 GeV) and m > \q±\ (heavy projectile quark, m c — 1.5 GeV). 
Radiation is suppressed for 9 < 9o due to the dead cone effect, whereby the reduction 
is stronger as the quark mass grows. The chosen parameters are pi — (0, 0, 10) GeV, 
p f± = (0.3,0.2) GeV, n = 0.5 GeV, and uj = 0.001 GeV, cj> = tt/2 for the emitted 
gluon. 

from the abelian total matrix rad^ CD = -^l o ? + ^?fo > i s shown as a function 
of the gluon emission angle 6, which encloses the directions of the gluon and the 
incident quark. For the same kinematics the cases of a light and a heavy projectile 
are contrasted and reveal significant suppression of the radiation pattern at angles 
9 < 9o, which is the dead cone effect. The heavy particle case is classified with the 
condition < m, and vice versa, hence it depends on the choice of the transverse 
momentum transfer q±. 

From analytical and approximated matrix elements one can derive the \/{uj 2 9 2 ) 
decrease of the radiation amplitude for 9 > 9q and a 9 2 /uj 2 slope for 9 < 9q, as can 
be seen in Fig.|2]for both cases. This behavior for 9 ~ 9q agrees with the dead cone 
factor. 

Some additional features seen in Fig. [3 are clarified in Fig. where the indi- 
vidual contributions for the light quark situation are depicted, including the non- 
abelian contribution. For angles 9 between 9 q = and 9 C = = ^ the 
pre-emission diagram dominates, and for 9 > 9 q pre- and post-emission interfere 
to yield a different color factor which accounts for the shift by | at 9 q in Fig. [5] 
The constant post-emission process for 9 < 9 C determines the constant radiation 
amplitude in this angular region. 

Although in Fig. [21 the picture of the radiation amplitude in the heavy quark 
case appears quite similar, here it is the pre- and post emission which interfere 
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Fig. 3. Contributions of individual diagrams to the total results in Fig. [3 for 
the light quark scenario. The three-gluon contribution is suppressed by orders of 
magnitude and so the full QCD result does not deviate from the abelian QCD 
radiation amplitude. 



with a magnitude of same order to obtain the dead cone behavior for 9 q < 9. The 
post-emission is dominant and constant for 9 < 9 q and thus the radiation amplitude 
becomes constant as well. 

Note, that in a corresponding QED situation the radiation probability will be 
further reduced for 9 > 9 q due to destructive interference, which in QCD is excluded 
by color factors. 

3.2. Non-Abelian Diagram and Factorization 

We continue to discuss to which degree the three-gluon vertex diagram lc con- 
tributes and how it may change the behavior of the abelian QCD single scattering 
radiation amplitude. In fact, the results for the kinematical situations considered 
so far for abelian diagrams are not effected by the additional three-gluon vertex di- 
agram. Figure |21 is not modified by this contribution; its smallness is demonstrated 
for the light quark case in Fig. |3] 

It can be shown that the non-abelian part becomes especially important for an- 
gles 9 > 9z g i u with #3gi u = ^%^, but this implies \q±_\ <u>, in order to consider such 
an angular configuration. This is for example realized in Fig. 0] where calculations 
of individual contributions in the heavy quark case are shown. It exhibits that the 
three-gluon contribution dominates in the whole 9 range. However, if one allows 
for 9s g i u < 1, this means to abandon the soft radiation limit. In other words, the 
condition |gj_| < u> implies that the momentum transfers, which are constrained 
by initial and final states, differ significantly in the radiation and elastic scattering 
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Fig. 4. Contributions of individual diagrams for a situation with a dominating 
three gluon vertex diagram. The parameters are chosen for #3 g i u < 1, so that 
03giu < is possible, m c = 1.5 GeV, p- = (0,0,10) GeV, p A = (0.03,0.02) GeV, 
/i = 0.05 GeV, uj = 1 GeV, <j> = tt/2. 

processes. But the latter serves as reference system for the definition of the radi- 
ation amplitude. The factorization, which is usually applied in the derivation of 
approximated matrix elements, becomes insufficient. In specific examples, this may 
be accounted for by an additional correction factor [GO]- Here we mention that this 
problem causes e.g. #o m Fig. to be shifted off the peak position. 

4. Validity of the Potential Model 

It is tempting to test the validity of the potential model beyond the assumptions 
of massless scattering particles. We analyse this while we contrast the scattering of 
an incident particle in the potential model with an analogue scattering of the same 
projectile on a target quark at rest. 

It was found that the radiation amplitude calculated in a quark-quark collision 
coincides well with the potential model result as long as the angles are small and 
the projectile is not heavier than the target. 

The results of a scalar QCD calculation for light targets and heavy projectiles 
are exhibited in Fig. [3] There the radiation amplitude in the potential model can 
be well-described by the application of a general dead cone factor, which is valid 
for arbitrary angles. For smaller 6 the potential model result becomes constant 
(cf. Fig. 2) and a simple dead cone factor does not longer apply. The dead cone 
suppression is weaker in a quark-quark scattering approach. This shows that for 
such an arrangement of colliding masses the potential model might be challenged, 
since the assumption, that no energy would be transferred to the target, is violated. 
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Fig. 5. The radiation amplitude |i?i| 2 as a function of the emission angle 6 is 
compared for the potential model, the quark-quark scattering on a light target and 
the general dead cone factor. The chosen parameters are rrid — 0.007 GeV (light 
quark), m c = 1.5 GeV (heavy quark), pi = (0,0,10) GeV, pi, Target = (0,0,0), 
q x = (0, 0.01) GeV, fi = 0, lo = 0.001 GeV, cf> = 3tt/2. 



This is intended to possibly point to the need of an improved model to describe the 
scattering processes in a deconfined medium for heavy projectile quarks. 

Besides this we have numerically justified the usual neglect of radiation con- 
tributions from the target lines in the A + -gauge. Also comparisons of scalar and 
spinor calculations did not reveal significant differences for small angles, light quarks 
and soft gluons. 

In the light quark scattering we find only configurations without spin flip to 
contribute to the radiation amplitude. In contrast, scattering a heavy on a light 
quark revealed that configurations become important where the spin of the light 
target is flipped. 

5. Conclusions 

This work is aimed investigating how the mass parameter of an on-shell quark 
passing a deconfined medium influences its radiative energy loss probability. A 
detailed discussion of the basic single scattering diagrams with one-gluon emission 
in the potential model was done. With the numerical calculation of the radiation 
amplitude in QCD we could disentangle interference effects, especially the influence 
of the non-abelian effects. 

The dead cone suppression factor was shown to emerge in specific angular re- 
gions, however the radiation amplitude deviates from this prediction for either small 
angles, due to the post-emission process, or in case of higher gluon energies lo. For 
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the latter possibility the three-gluon diagram, including the three-gluon vertex, be- 
comes responsible. We thus have confirmed the suppression effect due to the heavy 
quark mass, but find the dead cone factor (2) is not the correct modification in all 
kinematical situations. (This finding agrees with the one in [I14j. where an averaged 
dead cone suppression factor was also found to be insufficient.) We stress that not 
the mass parameter of the projectile itself but the ratio to the transverse momen- 
tum transfer is relevant to categorize light and heavy mass situations. Non-abelian 
effects become important for large gluon energies, that is to say greater than the 
transverse momentum transfer. 

We have indicated that the potential model might become inappropriate if 
heavy projectile quarks are present. In addition, it was emphasized that for consid- 
erably large values of the gluon energy one is confronted with factorization problems 
with respect to the elastic scattering part in matrix elements. In such situations 
the missing unique interpretation of the radiation amplitude questions whether it 
is sufficient to consider ratios of inelastic and elastic cross sections. 

The numerical confirmation of the LPM effect and the calculation of net energy 
loss of heavy quarks at finite energy along the lines of [ 1141 ITB*] are promising goals of 
further investigation, such as the extension of our approach to multiple scattering, 
i.e. higher orders in opacity, and multi-gluon emission. 
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